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Abstract. Full control over the spatio-temporal structure of quantum states 
of light is an important goal in quantum optics, to generate for instance single- 
mode quantum pulses or to encode information on multiple modes, enhancing 
channel capacities. Quantum light pulses feature an inherent, rich spectral 
broadband-mode structure. In recent years, exploring the use of integrated optics 
as well as source-engineering has led to a deep understanding of the pulse-mode 
structure of guided quantum states of light. In addition, several groups have 
started to investigate the manipulation of quantum states by means of single- 
photon frequency conversion. In this paper we explore new routes towards 
complete control of the inherent pulse-modes of ultrafast pulsed quantum states 
by employing specifically designed nonlinear waveguides with adapted dispersion 
properties. Starting from our recently proposed quantum pulse gate (QPG) we 
further generalize the concept of spatio-spcctral engineering for arbitrary X 
based quantum processes. We analyse the sum-frequency generation based QPG 
and introduce the difference-frequency generation based quantum pulse shaper 
(QPS). Together, these versatile and robust integrated optics devices allow for 
arbitrary manipulations of the pulse-mode structure of ultrafast pulsed quantum 
states. The QPG can be utilized to select an arbitrary pulse mode from a 
multimodc input state, whereas the QPS enables the generation of specific pulse 
modes from an input wavepacket with Gaussian-shaped spectrum. 
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1. Motivation 

Ultrafast pulsed quantum states of light play an increasingly important role in 
quantum information and quantum communication as they allow for efficient network 
synchronization and high data transmission rates. In general they feature a rich 
spectral mode structure, which is most naturally described in a broadband pulse-mode 
basis. This is not a new result in either classical or quantum optics [1]. For classical 
states all basis sets are formally equivalent and no specific choice can be distinguished. 
In contrast, it has been shown that pulsed quantum states of light exhibit an inherent 
pulse- mode structure, which is solely determined by their generation process [2]. 
Different kinds of applications require specifically tailored pulsed quantum states, be 
it single-mode states for linear optical quantum computation [3] or multimode states 
for high-capacity quantum information encoding. Thus, a thorough understanding of 
the spatio-spectral modal structure of ultrafast quantum states as well as the ability 
to exercise full control over that structure is an important goal in today's quantum 
optical research. 

In this paper we investigate the potential of engineered nonlinear waveguides 
for the manipulation of pulsed quantum states which cannot be achieved within the 
framework of linear optics. Special emphasis is put on an accurate description of 
the x^-process inside the guide, which takes into account rigorously the spatial 
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and spectral degrees of freedom. Thus quantitative measures can be derived for the 
efficiency of practical quantum optical devices. 

The paper is organized as follows. In section 2, we review the state-of-the- 
art methods of generating ultrafast pulsed quantum states and manipulating their 
inherent pulse-mode structure in bulk crystals and waveguides. We briefly discuss the 
latest developments and introduce new ideas by combining dispersion engineering 
techniques, which have become established by now for photon-pair preparation, 
with current methods of state manipulation utilizing x^ 2 ^-nonlinearities. In this 
context we analyse the experimental implementation of our recently proposed quantum 
pulse gates (QPG) [4] and extend the formalism further by presenting the concept 
of a quantum pulse shaper (QPS). In sections 3 and 4 we develop a theoretical 
framework for our devices. We start with the linear operator transformations for 
sum- and difference-frequency generation and derive the interaction Hamiltonian of 
these processes considering spatial and temporal degrees of freedom. Our analysis 
results in a completely quantitative model. Section 5 is dedicated to merging the 
derived theoretical framework with dispersion engineering methods known from state 
preparation, thus paving the way to real- world applications, the performance of which 
is investigated in section 6. Here we introduce realistic experimental parameters for 
our waveguide devices and demonstrate the capability to fully control the pulse-mode 
structure of ultrafast pulsed quantum states of light. Finally, in section 7 we highlight 
the most important results of this work, and end with an outlook on the use of QPG 
and QPS in continuous variable quantum information processing. 

2. Introduction 

In recent years different approaches have been introduced to prepare and manipulate 
ultrafast pulsed quantum states of light. One of the most common sources for 
the generation of photonics quantum states is parametric downconversion (PDC) 
in nonlinear crystals. This is mainly due to the rather simple experimental 
implementation of PDC sources and their ability to achieve high photon-pair 
generation rates. When pumped by ultrafast pulses, PDC processes generate pulsed 
bi-photons with broad spectra. However, these states are usually highly correlated 
due to the constraints imposed by energy and momentum conservation [5, 6]. Hence 
photon pairs are typically emitted into many inter- and intra-correlated spatial- 
spectral modes, the exact structure of which can be retrieved by applying a Schmidt 
decomposition to the biphoton amplitude distribution [7]. Upon detection of one of 
the photons, the other one is projected onto a mixed state of all possible modes, 
rendering it ill-suited for linear optical quantum computation applications [8]. The 
common way for overcoming this limitation has been narrowband spectral filtering to 
force the photons into one optical mode [9, 10]. However, this approach prohibitively 
lowers the photon generation rate as most of the generated signal is lost in the filtering 
process. It is thus not feasible for large scale quantum information applications [8] . In 
addition, only in the limit of infinitely narrow filtering one monochromatic, temporally 
de- localized mode is selected, and the photon's pulse characteristic is lost. 

2.1. Preparation of ultrafast pulsed quantum states with waveguides 

Only recently two new developments have made it possible to tackle the 
aforementioned problems. The use of integrated waveguide sources has a major 
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impact on the structure of PDC photon pair states. In a bulk crystal the generated 
photons are emitted at the natural phasematching angles. This poses two problems: 
Firstly, the collection of the pair photons is experimentally challenging and typically 
inefficient. Secondly, the pump field always couples to an infinite number of spectral- 
spatially correlated modes, and thus the probability to create a photon pair in one 
distinct mode becomes very low. In contrast, the emission in nonlinear waveguides 
is restricted to a well-defined set of discrete spatial modes defined by the waveguide, 
ideally allowing only the propagation of one individual mode in a single-mode wave- 
guide. \ It turns out that the probability of generating a photon pair in a distinct 
spatial mode is enhanced by several orders of magnitude [11, 12], since the total 
number of allowed modes is dramatically reduced inside the waveguide. Moreover it 
also leads to an effective decoupling of the spatial from the spectral degree of freedom, 
since any spatial-spectral correlation necessitates more than one spatial mode. Even if 
other spatial modes apart from the ground mode are guided in the waveguide, modal 
waveguide dispersion usually ensures that both photons in those modes are created 
at different frequencies. Thus they can readily be removed by applying broadband 
spectral filters on the output state [13]. 

The second step on the way to achieve complete control over the modal structure 
of the generated quantum states of light is spectral source-engineering. It has been 
proposed that, by choosing adapted dispersion properties, photon pair generation can 
be tailored such that signal and idler are emitted into one single spectral pulse mode 
each [8, 14]. Later this has been experimentally demonstrated for bulk crystal sources 
[15, 16, 17] and photonic crystal fibre sources [18, 19, 20]. Only recently this has been 
realized in a waveguided parametric downconversion source in a KTiOP04 crystal 
[21]. In this setup, the use of a waveguiding structure has led to an unprecedented 
brightness for sources of separable photon-pair states. Note that narrowband spectral 
filtering is not necessary with these sources as the generation process itself only allows 
one single spectral pulse mode. Thus, the generated photon pairs are genuine quantum 
pulses, and are completely separable, spectrally as well as spatially. 

2.2. Manipulation of the pulse-mode structure of ultrafast quantum states 

Up to now, the research on manipulation of pulsed quantum states has mostly 
been focussed on shifting their central frequency. It has been shown that sum- 
frequency generation (SFG) of single photons, in combination with subsequent 
photodetection, surpasses the efficiency of direct detection of near infrared single 
photons [22, 23, 24, 25]. Additionally SFG has been proven to conserve the quantum 
characteristics of the input photon [26, 27, 28] and it has already been utilized to 
implement measurement schemes with very high timing resolution, which overcomes 
long integration times of current single-photon detectors [16]. Only last year, SFG has 
been demonstrated for single-photon Fock states [29]. Note that recently also four- 
wave mixing in photonic crystal fibres has been employed to demonstrate coherent 
frequency translation of single photons [30]. This highlights the broad interest and 
the numerous application possibilities for these techniques. With more and more 
single-photon sources available in the visible range, difference-frequency generation 
(DFG) has now also attracted increasing interest. Recent experiments employ DFG to 

| Here wc do only consider guided modes neglecting any contributions, which could be present due to 
phase-matched substrate modes. This simplification is justified, because the continuous distribution 
of substrate modes can be easily filtered out by spectral or spatial filters. 
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Figure 1. Schematics of (a) a Quantum Pulse Gate (QPG) and (b) a Quantum 
Pulse Shaper (QPS). The QPG is based on engineered ultrafast sum- frequency 
generation in nonlinear optical waveguides. By shaping an ultrafast gating pulse 
one specific pulse mode from a pulsed multimode input state is selected and shifted 
to another frequency. Then it can easily be split off while leaving the rest of the 
state untouched. The QPS is based on engineered ultrafast difference-frequency 
generation in nonlinear optical waveguides. An input pulse mode can be converted 
into an arbitrary output pulse mode by nonlinear interaction with an ultrafast 
shaping pulse. The output pulse mode's shape is given by the mode of the shaping 
pulse. 



implement wavelength interfaces for quantum networks [31, 32, 33], which equivalently 
to the SFG process preserve the quantum characteristics of the input state. 

Despite this considerable progress, the generation of ultrafast pulsed quantum 
states with a specific pulse-mode structure, be it the number of excited modes or 
their shape, has not been explored, yet. Although the rich inherent mode structure of 
ultrafast optical quantum states is well-known, before the QPG [4] there has been no 
feasible way of controlling and manipulating the different modes separately. 

2.3. Quantum Pulse Gate and Quantum Pulse Shaper 

We now combine the findings from the field of quantum-state generation with 
the techniques from state manipulation. Applying source-engineering to frequency 
conversion reveals fascinating possibilities to achieve the desired goal of complete 
control over the pulse-mode structure of ultrafast quantum states. In [4] we have 
already proposed a quantum pulse gate (QPG), a device based on engineered ultrafast 
SFG in nonlinear waveguides. This device enables us to address different inherent pulse 
modes of an ultrafast pulsed quantum state of light individually as illustrated in figure 
1 (a). We would like to highlight that the QPG operation does not have any impact 
on the residual pulse-mode structure. This sets it apart from other experiments which 
focus on a direct manipulation of the spectral broadband-mode structure of ultrafast 
pulsed quantum states and employ pulse shaping of photon-pair states [34, 35, 36]. 
This alternative approach also leads to highly interesting results for entanglement 
based applications. Still, the manipulation is not pulse-mode sensitive in the sense 
of accessing and separating out a single-mode quantum state with specific temporal 
profile. In contrast, the QPG achieves mode selection by shaping an ultrafast, coherent 
gating pulse instead of the pulsed quantum state. The addressed mode is converted 
to the sum-frequency of input pulse and gating pulse and is thus easily accessible. 
In addition, different orthogonal pulse modes can be interconverted into each other, 
rendering interference between them possible. In this paper we elaborate on the 
QPG concept and come up with another fundamental device, the quantum pulse 
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shapcr (QPS) based on engineered ultrafast DFG. While the QPG addresses single 
pulse modes, the QPS enables us to convert an input quantum state with Gaussian- 
shaped spectrum into a single-mode quantum state with arbitrary shape (see figure 
1 (b)). Here, an arbitrarily chosen pulse form of the coherent shaping pulse defines 
the output pulse mode. We would like to mention that a similar idea of shaping 
quantum pulses by means of frequency conversion with dispersion matching has been 
proposed in [37]. In contrast to this earlier work we put special emphasis on the 
engineering of the dispersion characteristics of the used non-linear medium, such that 
single-mode operation can be ensured avoiding the insertion of any unwanted vacuum 
contributions. 

Using QPG and QPS, pulsed quantum states can be generated and selected 
with unit efficiency in arbitrary pulse forms and encoding of quantum information 
in broadband mode basis and the successive read-out become possible. Therefore 
QPG and QPS will enable the implementation of quantum communication protocols, 
which exploit the rich pulse-mode structure of ultrafast states. 

3. Linear transformations for SFG and DFG in comparison with PDC 

In this section we qualitatively discuss the nonlinear three-wave mixing processes 
SFG, DFG and PDC, highlighting their formal similarities as well as examining 
their differences. In such a three-wave mixing process, three electrical fields interact 
inside a nonlinear medium, and the interaction Hamiltonian in the rotating-wavc 
approximation is of the form 



The E\ ' (r, t) describe the positive frequency parts of the interacting electric fields, 
X^ 2 ' is the second order nonlinearity of the medium. In PDC and single-photon SFG 
and DFG, two of the three fields are generally considered quantum mechanically. The 
remaining field is a bright, immutable pump field which can be treated classically. In 
this case the interaction Hamiltonian becomes bilinear and Heisenberg's equation of 
motion yields linear input-output transformations for the creation and annihilation 
operators. Depending on which of the three fields is defined as pump, one can 
distinguish two flavours of processes which are characterized by different linear 
operator transformations. This can be derived when considering a single-mode 
approximation to equation 1: 



Firstly we assume that field E a (r, t) is the classical coherent pump field. We insert its 
classical amplitude a into the above equation and find 



The resulting operator formally corresponds to a two-mode squeezing operator 
(compare e.g. [38]), which means that this case describes PDC. Depending on the 
pump power of the bright field, either the photon pair characteristics (low power 
regime) or the squeezer characteristics (high power regime) dominate the PDC output 
state. The linear transformations between input and output operators evaluate to 




(1) 



Hi nt oc atf & + a' be. 



(2) 



i?int oc a trc' + a* be. 



(3) 



b^ cosh(C)6 - sinh(C)c t , 
c — s- — sinh(£)&T + cosh(C)c, 



(4) 
(5) 
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(a) (b) (c) 




Figure 2. Schematics diagrams illustrating two different flavours of second-order 
nonlinear processes. All processes are pumped by a classical undcplctcd field at a 
frequency lo v , the dashed arrows indicate vacuum modes, (a) In a PDC process 
two photons get created and the evolution operator for this process is a two-mode 
squeezing operator, (b), (c) In a SFG or DFG process one incoming photon is 
annihilated and and outgoing photon at another frequency is created. Assuming 
that the frequencies involved in the two processes are equal, one can readily see 
that SFG and DFG are similar, yet reversed processes which is indicated in the 
schematics by the different directions the arrows point to. The corresponding 
evolution operator for these processes is equivalent to a beamsplitter. For further 
details see text. 



where the parameter £ depends on the pump power and is related to the amount of 
squeezing in the generated pair state. This is discussed, for instance, in [39]. 

We find the other flavour of processes by assuming that field E\,(f,t) 
corresponds to the pump field. We substitute its classical amplitude /? in (2) and 
obtain 

H mt oc /3 act + p* ate. (6) 

This expression is formally equivalent to an optical beamsplitter Hamiltonian and we 
can use the well-known beamsplitter input/output transformations for the operators 
a and c. 

a — > cos(0)a — isin(0)c, (7) 
c ->■ — isin(0)a + cos(#)c. (8) 

We identify 9 with the beamsplitter angle which depends on the pump power 
and the strength of the nonlinear interaction. This will be discussed later in great 
detail. We interpret this x^-process as a beamsplitter which diverts optical beams 
into different frequency output ports depending on their initial frequency. Note that 
in single-photon quantum optics this Hamiltonian describes SFG as well as DFG. In 
classical nonlinear optics, however, DFG is understood as a stimulated process. The 
bright pump field has the highest frequency and the process is seeded with a weak 
input field which gets enhanced through continuous conversion of pump photons. In 
that case similar operator transformations as for PDC are valid and the process could 
also be interpreted as seeded PDC. In contrast, we assume a single (or few) photon 
input state, which has the highest frequency and the "seed" field is the bright field. 
By pinning the pump field to a fixed value in our process model, we exclude that 
stimulation can occur and the process becomes formally equivalent to SFG. Note that 
the usual no-pump- depletion approximation (dE v /dz = 0) needs to be interpreted as 
a no-pump- enhancement approximation in this case. These findings are schematically 
depicted in figure 2. We move on to the derivation of the interaction Hamiltonian for 
quantum mechanical frequency conversion inside an optical waveguide. 
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4. Quantitative derivation of the SFG and DFG interaction Hamiltonian 

4-1- Spatial mode considerations in a monochromatic approach 

The interaction Hamiltonian of a frequency conversion process can be expressed as 

H int = -d eS e J d 3 rE p (r,t)El + \r,t)E(-\r,t) + h.c, (9) 

where d c ff denotes the effective nonlincarity, E p is the classical pump field and e[ +) 
and Eo ' denote the operator expressions for the input signal and the converted 
output, respectively. As the interaction happens inside a nonlinear optical waveguide, 
the propagation of the fields is restricted to one direction, which is given by the 
waveguide axis and which we define as z-direction. The ultrafast pump field then 
reads 

E p (r,t) = A p f p (x,y) [ d^a^c"^^. (10) 



Here, a(w p ) is the normalized spectral amplitude of the pump. The function f p (x,y) 
describes the transverse spatial distribution of the pump field with J d 2 r|/ p (a;, y)\ 2 = 1 
and /3 P is the propagation constant of the corresponding transverse mode. By requiring 
that the area integration over the field intensity I = ^cn p eo\E\ 2 , where n p denotes 
the refractive index at the pump frequency, corresponds to a power, we find that the 
amplitude A p is related to the average pump-pulse peak power P p by 

( 2P p \ 1/2 

Ap \ce Q n p (u p )\ J doj p a(u! p )\ 2 J 

We implicitly make use of the slowly-varying envelope approximation in this 
calculations, which is valid as we consider only pulses with Aui <C loq. This also 
means that we can neglect the frequency dependence of the propagation constant /3 P 
in (10). 

To derive expressions for the quantized fields in a nonlinear waveguide we start 
from the electric field operator for a propagating field in a dielectric with finite cross- 
section area A, given in [40] . Note that we assume the wavevector components k x and 
k y of the quantum field to have fixed, finite values. 

1 /2 

E^(x,y,z,t)=i fdu( —, r j ] a(uj)e- iujt+i( - k * x+k yy +k * z l (12) 

J \A7re cn{u})AJ 

In a nonlinear waveguide with field propagation along z-direction, the solution of 
the Helmholtz equation yields a discrete spectrum of valid propagation constants 
Pmn = kz and a set of allowed, localized transverse modes {f m n{x, y)}, determined 
by the boundary conditions of the guiding geometry. The indices m and n denote 
the order of the transverse mode in x- and y-dircction. Each (5 mn corresponds to 
exactly one mode and, in the case of a symmetric situation, the /3 mn for corresponding 
modes (e.g. /3oi and /3io) become degenerate. For better readability we only use 
one index m to describe the modes. The spatial localization of the fields implies 
continuous spectra of the individual wavevector components k^ and k y m \ given 
by fm(k x ,k y ) = TT [fm(x, t/)]§. The electric field inside a nonlinear waveguide is 

§ We deploy the symmetric definition of the Fourier transform that is f(k) = ^ f dxf(x)e~ lkx 
and accordingly f(x) = * J dkf(k)e lkx . 
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accordingly comprised of a superposition of quantum fields of the form given in (12), 
each corresponding to a single propagation direction. Summing over the transverse 
waveguide modes and integrating over the k^ and k y we find: 

E(+\x,y,z,t)=iJ2 I ( HU y^He^^ x 
^ J \&ir£ocn m (u) J 

x [dk^Uk^f m (k x ,k y )c ik * m)x+ik( y m) y (13) 



This can - due to the Fourier relationship between position and momentum - be 
written as 

E^(x,y,z,t)=iJ2^fm(x,y) f du (- ^- r ^] 1 a m Mc-- t+i ^. (14) 

J \AirE cn m {u) J 

We would like to point out that we account for the cross-section area by the spatial 
distributions f m (x, y) which are normalized such that J dxdy\f m (x, y)\ 2 = 1 and which 
have units of inverse meters. Moreover we assume that within the frequency range 
of the considered fields the variation of the spatial properties is negligible, due to the 
narrowband approximation Acj <C u>q. We substitute the electric field operators into 
(9) and rephrase the interaction Hamiltonian for single photon frequency conversion 
as 

f d c ghn ^ / WiW 
H ln t = A p * / : 



v *H,i'*'o,m 



x / dxdyf p (x,y)fu(x iy )f* m {x,y) x 
x J tejV,±b,*Po, m )' x 

x /d Wp d Wi d Wo a( Wp )e- i ^ ±w '^°) t a ; ( Wi )4( Wo )+h.c, (15) 



where we discriminate between SFG and DFG. Here labels i and o denote input and 
output fields, whereas the indices I and m describe the transverse spatial modes of 
input and output field, respectively. 

Now we move on to the calculation of the time evolution of the input quantum 
state during the conversion process. Note that the interaction of the Hamilton operator 
of equation (15) is time-dependent and thus the exact solution has to take into account 
time-ordering effects. Here we present an approximate solution which neglects time- 
ordering effects, in order to emphasize the conceptual structure and to illustrate the 
main idea. In 7 we validate this approach by comparing the approximate solution with 
rigorous calculations we performed. We find that the shape of the mode functions does 
not change significantly when taking into account time-ordering, but the maximum 
conversion efficiency drops to 90%. Still, these findings confirm that the analytical 
solution leads to reasonable results and can safely be applied. Hence we write the 
time evolution of the quantum state during the conversion process 

Wont = U{t) |V) = exp (-i y dtH int (t)j |V')o ■ (16) 

Thus, we need to perform the time integration of the interaction Hamiltonian given 
in (15). This is a well-known procedure discussed for PDC in great detail in [8]. We 
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only present the result here as the calculation, including the waveguide aspects, is 
straightforward . 

2d cS hTT 2 



dtH int (t) 



-A P LJ2 



UJiLU 



1 



x / dwidw a(ctii )^ iTO (cJi, (j„)ai(wi)4K) + h.c. 



(17) 

Here, L is the length of the nonlinear waveguide. The function a(u>i ) is the spectral 
pump distribution defined as a(uj — w;) for SFG and a(ojj — lo q ) for DFG, respectively, 
whereas the function (j)i.m {^U w o) characterizes the phasematching and is given by 



0Z,m(Wi, Wo) 



sine 



exp 



-0.193- 



A/3, m L 



(18) 



The expression A/3; m describes the phase-mismatch of the propagation constants and 
evaluates to A/3 z , m = /3 p + /3 M - /3 0)TO - — for SFG and A A,m =P P ~ A,/ + - x 
for DFG, respectively. Finally, A is an optional poling period for quasi- phasematching 
inside the waveguide. Following the usual conventions we define an effective interaction 



A(eff). 



dxdyf p (x,y)fi t i(x,y)f* (x,y) 



(19) 



Note that this should not be mistaken as a geometric area defined for instance by 
the waveguide cross section. Instead it describes the overlap of the transverse spatial 
modes of the three interacting fields inside the nonlinear waveguide. This result also 
implies that simply using a smaller waveguide - while not changing the modal overlap 
characteristics - will not alter A^- e ^^ and will therefore not have any impact on the 
efficiencies of the processes. The product of pump distribution and phasematching 
function is conveniently defined as joint spectral distribution function 

G; !TO (wi, ui ) = ——a(u;i )<f)i <m (u}i,cj ), (20) 

which describes the mapping between input and output frequencies for a 
specific pair of spatial modes l,m. The normalization factor N^ m reads 

(/ dwidw o |a( io )0 ;im (wi, w )| 2 ) 1/2 . 



4-2. Broadband pulse mode picture 

The description derived so far has been in terms of monochromatic creation and 
annihilation operators. However, since we concentrate on x^-iuteractions between 
ultrafast pulses, a much more natural approach is to consider broadband pulse modes. 
A suitable pulse-mode basis is found by applying a Schmidt decomposition to the joint 
spectral distribution function: 

G l>m ( WijWo ) = £>f m V?' m W?' m Vo)- (21) 

3 

Equation (21) yields two correlated sets of orthonormal broadband pulse- mode 
functions {<^' m )(wi)} and {-0(^ m ) (lu )}. The diagonal values K^f' m ^ are the real and 
positive Schmidt coefficients and satisfy J2j( K< j' m ^) 2 = 1- ^ i s wcu known for PDC 
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that the basis sets of the Schmidt decomposition, and thus the modal structure of the 
photons, are uniquely defined [2] . The same argument can also be applied here in the 
context of SFG and DFG. As for PDC we define broadband creation and annihilation 
operators: 

Aj,i, m = duJiipf' m \uJi)ai(uji), (22) 



Cj,i, m = J duj ^ l ' m \iLi )c m (uj ). (23) 

Substituting those, we rewrite the expression for the time-integrated interaction 
Hamiltonian from (17), and arrive at the broadband pulse-mode picture 

2d cS tm 2 ^ / WiW N L . 



, V Th,Wo,m / A<- e ff) 



3 

6 ^m (-Ij.n.r.,,,, + A] lm C^ m ) ■ (24) 



By introducing the effective coupling constant 9 3t i tW , in (24), we reveal the simple 
beamsplitter structure of the Hamiltonian [41], as already announced in (6). In 
contrast to a conventional beamsplitter however, this Hamiltonian does not couple two 
fc-modes (or beam paths), but rather two broadband pulse modes Aj t i >m and C^i_ rn 
at different frequencies! This is a unique feature of ultrafast frequency conversion 
processes and makes them ideal candidates for the implementation of the QPG and 
QPS. 

5. Pushing towards applications 

5.1. General non- engineered SFG and DFG 

QPG and QPS are unique in their single-mode operation on broadband pulse modes. 
In this section we discuss the implementation of genuine QPG or QPS in a feasible 
experimental setup. We restrict the analysis to only one pair of transverse spatial 
modes (I, m), which simplifies the notation but does not change the underlying physics. 
In the experimental setting the selection of one spatial mode can be accomplished by 
broadband spectral filtering [13]. In this case the time-integrated, effective SFG- and 
DFG-Hamiltonian from (24) reads 

f dtH int (t) = hJ2 0, (A 3 C] + A]Cj) , (25) 

with the broadband operators defined as 

Aj = f du)i<pj(u)i)a(u)i), (26) 

Cj = I dio ipj{io )c(uj ). (27) 



The complete, orthonormal function sets {(pj(uj[)} and {tpj(ui )} represent the intrinsic 
pulse-mode structure of the SFG- or DFG-process, obtained from the Schmidt- 
decomposition of the joint spectral distribution function G(u>i,u> ) = a(LJi )<f>(wi, ui ). 
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Figure 3. Conversion efficiencies n--, joint spectral distribution functions 
G(b}i,L0 o ) and Schmidt coefficients Kj for the first four pairs of pulse modes 
(J = 0...3). (a) Non-engineered process with several Kj 7^ 0. For a given 
overall beamsplitter angle 9, all modes are converted to a certain extent given 
by r]j = sin(Kj ■ 9) 2 . However, an overall unit efficiency can generally not be 
accomplished, (b) Source-engineered process with one predominant Kj fa 1. By 
choosing an appropriate 8, pulse mode <p(uJi) can be converted into pulse mode 
ip(uJo) with unit efficiency, allowing for QPG operation. 



On the one hand these are determined by the pump pulse characteristics, but on 
the other hand they also critically depend on the nonlinear waveguide's material and 
modal dispersion properties. We have already stressed the formal equivalence between 
the expression from (25) and a sum of optical beamsplitter Hamiltonians. Hence the 
linear transformation for the broadband operators can be readily written as 

.1, >ens^;.i ; isin^,)^, (28) 

corresponding to a pulse mode conversion between tpj(uj\) and ipj(co ) with efficiency 
rjj = sin 2 ((9j). According to (24) the coupling constant 9j is given by 

2d cS ir 2 LN [2 2^ / P p ~ _ 

'~ Ky c \jce n p n- ino \SAu p a{u v )\^ AW) ~ "fO-W) 

Here, 9 is an overall beamsplitter angle defined by the process parameters. Its impact 
on the different modes j is given by 9j, where, for each mode, the overall beamsplitter 
angle is weighted with the corresponding Schmidt coefficient Kj. In figure 3 (a) we 
illustrate a general, non-cnginccrcd SFG. We show the joint spectral distribution 
function G(uii,u) ) as well as the Schmidt coefficients Kj for the first four pairs of 
pulse modes and plot the conversion efficiencies r\j versus the beamsplitter angle 9. It 
is obvious that, for any given value of 9, all pulse modes with Kj ^ are converted 
to a certain extent. Yet, in general neither single-mode operation is achievable, nor 
can conversion with rjj = 1 for different pulse modes simultaneously be accomplished. 
We note an exception to this rule: Under certain conditions (e.g. a cw pump) input 
and output modes are perfectly correlated. Then all Kj share the same value and all 
modes are converted with the same efficiency. The process is then highly multimodc 
but the overall efficiency can reach unity for high pump powers. 

5.2. Source- engineered SFG and DFG - towards genuine QPG and QPS 

We have shown that SFG and DFG in general are multimode processes. But for 
QPG and QPS we require single-mode operation in order to avoid signal degradation 
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introduced by vacuum contributions and to achieve unit efficiency. Reducing the 
intrinsic pulse-mode structure of a x^-nonlinear process to only one pair of modes has 
been extensively studied in PDC, where source-engineering led to the desired results 
[8, 15, 16]. Experimentally this is accomplished by group- velocity matching inside the 
nonlinear medium. If pump and either signal or idler share the same group velocity, 
the phascmatching function becomes parallel to one of the axes when plotted in an 
(w s , oji)-diagram. Then, the Schmidt decomposition yields - given that the process is 
pumped by an ultrafast pump - only one pair of pulse modes which is excited with 
unit efficiency. 

We transfer this insight to our analysis of SFG and DFG and employ it for 
spectral engineering of the conversion. The time-integrated Hamiltonian from (25) for 
our special case reduces to 

J dtH iDt = M(i(7 t + i f C*), (30) 

and can be interpreted as a beamsplitter operating on only one pair of pulse modes 
<p(u>i) and ip(co ). As an example, we show an engineered case in figure 3 (b) where 
pump and input signal are group-velocity matched. Note that, in contrast to the 
previous non-engineered case, the joint spectral distribution function now shows no 
spectral correlations between input and output frequencies that is it is oriented along 
the axes of the diagram. As in PDC this is a direct consequence of the horizontally 
oriented phasematching function and thus of the group-velocity matching. The 
distribution of the Kj reveals that only one coefficient kq differs significantly from 
zero. This is also reflected in the plot of the conversion efficiencies. Only one single 
pulse mode is addressed and, by choosing 9 = converted with unit efficiency and 
no vacuum is coupled into the signal beam. Hence we find that, by group- velocity 
matching pump and either input or output, we can achieve genuine single- mode 
operation and therefore implement QPG and QPS. Note that the data in figure 3 
have been calculated using our modeling and realistic experimental parameters, which 
are specified in section 6. 

Knowing how to achieve single- mode operation of SFG and DFG, the next step 
is to investigate how we can exact complete control over the pulse modes y>(u)i) and 
ip(u) )- A QPG selects a specific pulse mode from an input state and a QPS generates 
an arbitrary pulse mode from a Gaussian input mode. Hence, for QPG we require 
control over </?(wi), whereas for QPS we require shaping of ip(ui ), respectively. In 
figure 4 (a)-(c) we consider QPG. Shown are the phascmatching and pump functions 
as well as the resulting joint spectral distribution function. Note that the axes are 
given wavelength units rather than frequency for convenience. We find that the output 
mode tp(u) ) is defined solely by the phasematching. We performed calculations for 
three different spectral shapes of the pump and it is obvious that the input mode 
<p(u)i) has the form of the respective pump mode. Thus, in a QPG, spectrally shaping 
the bright gating pulse leads to the selection of an arbitrary pulse mode. In contrast, 
figure 4 (d)-(f) illustrate the situation for QPS. Now, pump and output are group- 
velocity matched, causing a vertical phasematching function. Again we calculated 
three different spectral pump distributions. We find that the spectral shape of the 
pump now defines the output mode ip(uj ) of the QPS and the input mode <p(ui) is 
fixed by the phascmatching. Summing up these findings, we end up with the following 
correspondences: 

a(w p ) -t <p(tJJi), (f>(oJi, uj ) ->• i>(oj ) for QPG, (31) 
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Figure 4. Joint spectral distribution, pump and phasematching function as well 
as pulse modes <*p(u\) and V(aj ) for QPG (a)-(c) and QPS (d)-(f). Note that 
the functions are plotted against wavelengths for reasons of convenience. In a 
QPG, the intrinsic pulse mode ip(u>i) can be manipulated by shaping the bright 
gating pulse, whereas pulse mode ^(lu ) is fixed by the phasematching function. 
Therefore, arbitrary input modes are mapped to the same output mode, allowing 
for interference of formerly orthogonal modes. Contrary to that, in a QPS, shaping 
the bright pulse defines the output pulse mode \j)(u>o). The input pulse mode is 
now defined by the phasematching function. Thus, an arbitrary mode can be 
generated from an input which is matched to ip{u)\). The data presented here is 
calculated using realistic experimental parameters, specified in section 6. 



<f>(ui,u ) ->• (p(oJi), a(w p ) -> tp(uj ) for QPS. (32) 

We demonstrated that we can achieve complete control over the required pulse mode 
of the QPG or QPS by shaping the bright gating pulse or shaping pulse, respectively. 
Note however that we considered only the intrinsic modes of the QPG and QPS, which 
do not necessarily have to coincide with the pulse-mode structure of an input signal. 

5.3. Mode matching a QPG or a QPS 

Given a specific QPG or QPS, the input signal's pulse-mode structure must coincide 
with the pulse-modes {(p(wi)} accepted by the device, in order to guarantee mode 
selectivity and high conversion efficiency. We first discuss this for the QPG. We 
have shown that the QPG can be easily adapted to a wide range of input signals by 
spectrally shaping the coherent gating pulse. The output mode V>( w o) is solely defined 
by the phasematching function and is independent from the pump pulse shape. It can 
typically be approximated by a Gaussian spectrum[8]. Hence any selected mode from 
an input state is mapped to the same output mode. 
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Figure 5. (a) Mode-matched QPG. The bright gating pulse has the same 
duration as the input signal, leading to the matching of the input pulse-mode 
structure and f(oJi) of the QPG. Only a single pulse mode of the input - the 
one which overlaps with <p(uii) - is selected and converted with unit efficiency, 
(b) Mode-mismatch in a QPG. The gating pulse and input signal have different 
durations, leading to an overlap of ip(u>i) with all input signal modes with the 
same parity. Hence, all of those modes are selected and converted to a certain 
degree. QPG-opcration is then not possible. 



In figure 5 (a) we illustrate this situation. Input signal and gating pulse share 
the same duration and the QPG is mode-matched to the input. Only the desired 
mode from the input signal is selected and converted with unit efficiency. In contrast, 
figure 5 (b) demonstrates the impact of a mode mismatch on the QPG operation. 
The gating pulse duration significantly differs from the input pulse duration and the 
intrinsic QPG pulse mode f{uJi) overlaps with all signal modes of the same parity. 
We end up with a case similar to multimodc SFG, with the only difference that the 
diverse conversion efficiencies for the modes are due to the different overlaps between 
and the corresponding signal-state modes. Because an overall efficiency of unity 
can not be achieved here and the process is not mode-selective anymore, it becomes 
clear that careful mode-machting is vital for a successful QPG implementation. 

The situation is different for QPS: The phasematching function is vertical in the 
(cl>;, cj )-plane (compare figure 4 (c)-(f)) and the input mode ip(uji) is now defined 
by the phasematching function alone. The QPS accepts only Gaussian input modes 
which are matched to However, shaping the bright pulse allows for defining 

the output mode ^(w ), rendering it possible to generate any pulse mode from an 
input pulse with a Gaussian spectrum. If the input state is not matched to this 
does not change the spectral form of the output pulse. The drawback is that not the 
complete input gets converted and vacuum contributions are introduced. 

We note that the same physical nonlinear waveguide device could be used as QPG 
or QPS, depending on whether SFG or DFG is implemented: this shows that QPG 
and QPG can be seen as reverse operations of each other. The results illustrate that 
QPG and QPS are versatile tools which can be easily adapted to a large range of input 
and output states, making them highly flexible and appealing for many applications. 



6. Performance of QPG and QPS considering realistic experimental 
parameters 

We conclude our analysis demonstrating the experimental feasibility of a QPG and 
derive, with the help of the theoretical model outlined in sections 3 to 5, an expression 
for the pump power for maximally efficient operation. The results of the calculation 
apply to QPS as well, since both devices can be implemented in the same nonlinear 
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Figure 6. Group velocities of the fundamental waveguide modes in a Ti- 
indiffused PPLN waveguide. For a 1550nm input oriented along the ordinary 
axis, an extraordinarily polarized pump pulse centered around 870nm is group- 
velocity matched. For a QPG, this consequently leads to an output at 557nm 
which has to be oriented along the ordinary axis. It becomes obvious that for a 
wide range of input signal wavelengths, group- velocity matched gating pulses can 
be found which still satisfy feasible experimental parameters. 
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Figure 7. Transverse spatial mode profiles of the fundamental waveguide modes 
in a Ti-indiffuscd PPLN waveguide. The mode profiles were calculated using 
a finite elements method. The overlap between the pump and the output 
exceeds 99%, because the guiding for ordinarily polarized fields is not as strongly 
pronounced as for extraordinarily polarized fields. 

waveguide. SFG phascmatching implies an existing DFG phasematching, only the 
roles of input and output field are interchanged. The coupling constant 9 is the same 
for both processes. The bright pulse, used as gating or shaping pulse depending on 
the application, will be called pump in this paragraph for ease of reading. The input 
field is at 1550nm and the third field is referred to as output. 

As a key point for the experimental setup we require that it can be operated at 
1550nm. The constraint of group- velocity matching determines the pump wavelength 
as soon as the input wavelength gets fixed, which in turn then also defines the output 
wavelengths due to energy conservation. We assume that the conversion takes place 
in a Ti-indiffuscd PPLN waveguide with a length of L = 10mm and at a temperature 
of T = 190°C to prevent the impact of photorefraction. Effective Sellmeier equations 
for the three participating fields were obtained by calculating the effective refractive 
indices of ordinary and extraordinary polarized fields with a finite-element method 
described in [42]. The effective equations were then fitted against the calculated 
values. Note that the following calculations are based on these effective Sellmeier 
equations. 
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In figure 6 we plot the group velocities for the ordinary and extraordinary crystal 
axes and assume for our modeling that all fields propagate in the fundamental 
transverse waveguide mode. If the input light is ordinarily (TE-) polarized with 
central wavelength of 1550nm, we find that the group- velocity matched pump has to be 
extraordinarily (TM-) polarized and centered around 870nm. The ordinarily-polarized 
output is then at 557nm. From figure 6 we can clearly recognize that a group- velocity 
matched pump can be found for any input, as long as the input is o-polarized. The 
effective refractive indices of the participating fields calculate to n p = 2.18, m = 2.21 
and n = 2.32 and we derive a periodic poling period of A w 4.28/im required for 
quasi-phasematching inside the waveguide. 

In figure 7 we plot the transverse spatial distributions of the input, pump and 
output modes, also obtained with the finite-element method from [42]. From these 
we calculate the effective interaction area « 64^m 2 . The conversion efficiency 

for a single- mode operation is r] = sin 2 (#) and the condition for unit efficiency can be 
specified by 



Assuming an input pulse duration of roughly 300fs, we calculate a required pump 
peak power of P p « 22W for optimal conversion efficiency. If a pump laser system 
with a repetition frequency of 76MHz is used, we obtain an average pump power of 
P av w 0.5mW inside the waveguide. This leads to required average pump powers 
of a few mW in front of the QPG or QPS, taking into account realistic waveguide 
coupling losses. These values are lower than formerly reported pump powers for 
similar experiments [22, 23, 24, 27, 29], owing to the careful source-engineering we 
applied to our process. This grants a significant advantage over experiments without 
spectrally engineered SFG, even though we employ a cross-polarized process with an 
effective nonlincarity that is lower by an order of magnitude compared to a process 
where all three fields arc oriented along the extraordinary crystal axis. 

7. Conclusion and Outlook 

In conclusion, we presented a feasible way to achieve complete control over the pulse- 
mode structure of ultrafast pulsed quantum states of light. We combined findings 
from quantum state generation and techniques from state manipulation, by applying 
spectral-source engineering and integrated optics to frequency conversion of ultrafast 
single photons. We showed that single-mode ultrafast sum- and difference-frequency 
generation in x^'-nonlinear materials are possible and analyzed two highly flexible and 
versatile devices, namely quantum pulse gate and quantum pulse shaper. The QPG 
is based on ultrafast SFG and offers the possibility to select arbitrary pulse modes 
from an ultrafast multimode input state. The selected mode gets converted with unit 
efficiency and is mapped onto a Gaussian output mode. The residual mode structure 
of the input is left intact, allowing for cascaded operation to convert multiple modes. 
As all input modes are mapped onto the same output mode, interference of formerly 
orthogonal states becomes possible. In contrast, the newly-introduced QPS is based 
on DFG and implements the reverse operation of a QPG. It enables us to generate an 
arbitrary pulse form from a Gaussian input mode. The output mode is defined by a 
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bright shaping pulse, thus highly flexible state preparation can be achieved. We have 
presented a quantitative analysis of QPG and QPS and derived feasible experimental 
parameters with which the proposed devices can be implemented, rendering them 
practical instead of merely conceptual. 

As a final remark we would like to point out that our analysis is in no way 
constrained to single-photon states. Although we consider single-photon input states, 
the introduced concepts can be generalized to classical and non-classical multi-photon 
states. In this framework the use of QPG and QPS provides the attractive opportunity 
to successively select and spatially separate arbitrary pulse modes from a multimodc 
input state while leaving the residual beam intact. Employing a series of QPGs 
operating on the same pulse mode in each arm of a multimode twin-beam squeezer 
source allows for a feasible implementation of non-Gaussian operations and thus 
constitutes an important step towards the realization of multimode continuous variable 
entanglement distillation. QPS on the other hand can be used to synthesize multimodc 
continuous variable Gaussian states in a mode-by-modc fashion. The prepared states 
can then subsequently be transmitted as a bundle, since they do not interact with each 
other and they all experience the same dispersion during transmission and therefore 
stay orthogonal. This paves the way towards dense channel-multiplexing in continuous 
variable quantum communication. 
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Appendix. Impact of time-ordering 

Since the interaction of the Hamilton operator of equation (15) is time-dependent, it 
might be assumed that time-ordering effects have a major impact on the intrinsic mode 
structure of the process, in particular if a perturbative solution is not sufficient. This 
case is associated with unit conversion efficiency, needed for perfect QPG and QPS 
operation. In our analysis we have nevertheless deployed the approximate solution 
which neglects time-ordering. The impact of time-ordering on the process of ultrafast 
PDC has been thoroughly investigated in [43], with the result that time-ordering 
mostly affects the amplitudes but not the shapes of the intrinsic pulse modes. In [37], 
the authors actually study a three-wave mixing process and find in their numerical 
simulation no major discrepancy with their analytical solution. This already indicates 
that, at least for low conversion efficiencies where a perturbative solution is sufficient, 
time-ordering can be neglected. However, since we aim for conversion efficiencies of 
unity and cannot conclude for sure that the above results remain valid in our case, we 
performed rigorous numerical simulations which take into account all time-ordering 
effects. Note that the results presented are all obtained for the case of maximum 
conversion efficiency. Additionally, the simulated processes are the ones discussed in 
this work. That is, pump pulse and input pulse are group-velocity matched and the 
processes have dccorrclatcd joint spectral distribution functions as is the case in figure 
3(b). 
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Figure 8. Input modes ip(u)i) as well as output modes ip(u) ) of the considered 
processes, obtained with the analytical solution (blue) and the rigorous calculation 
(red), respectively. Obviously, time-ordering has only a small impact on the actual 
mode shape. This can simply be corrected for in our schemes by adjusting the 
spectrum of the bright pump pulse. Note that the oscillations in the output modes 
originate in the sine function which describes the phasematching. 
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Figure 9. Schmidt coefficients obtained from the analytical as well as the 
rigorous solution. Time-ordering causes a slight shift in the excitation of the 
different modes. In addition it leads to a drop in the maximum conversion 
efficiency to 90%, as compared to a unit conversion efficiency reached in the 
analytical solution. 

Figure 8 depicts the analytical as well as the rigorous input modes y a na(wi) and 
</5rig(wi) and the corresponding output modes V'ana(w ) and ipn g (uio), respectively. It 
nicely illustrates that time-ordering has only a slight impact on the shape of the modes, 
as expected from [43] . The change in the modes can easily be compensated for in our 
proposed scheme, by adjusting the spectrum of the bright gating pulse. Note that 
the oscillations in the output modes originate from the sine function which describes 
the phasematching. These also cause the slight multi-modeness which can be seen in 
the Schmidt coefficients in figure 9, where the first higher order mode is also excited 
with a certain probability. Comparing again analytical and rigorous solutions, we find 
that time-ordering slightly shifts the weights between the different modes. However, 
we want to note that no new modes occur in the process due to time-ordering. The 
main difference is found when considering the maximum conversion efficiency. It turns 
out that this value drops in the rigorous solution to 90% instead of the unit efficiency 
obtained with the analytical approach. 
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This behavior will be thoroughly analyzed in [44]. We want to stress here that 
all characteristics introduced by the side lobes of the sine function can be washed 
out through careful design of the nonlinearity inside the waveguide as shown in [45] . 
Therefore the device performance calculated here only represents a lower bound and 
might be increased in the future, for instance by implementing a Gaussian shaped 
phascmatching function. 

In conclusion we find by comparing analytical and numerical solutions that 
the assumption that time-ordering can be neglected is, in fact, a rather good 
approximation, even for the cases of high conversion efficiencies analyzed in this work. 
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